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SUMMARY 


Linearized theory for compressible unsteady flow is used to treat 
the problem of a partial-span rectangular control surface, hinged at its 
leading edge and mounted on a rectangular wing, oscillating in supersonic 
flow. The motions of the wing-aileron configuration are assumed to consist 
of vertical translation, pitching, and aileron rotation. The velocity 
potentials for this conf iguration are derived as power series in terms 
of a frequency parameter to the fifth power of the frequency of oscilla- 
tion. From these potentials closed expressions for section force and 
moment coefficients at any spanwise location are derived. The section 
force and moment coefficients are employed in a particular flutter study 
in which the spanwise position of the aileron is varied. 


INTRODUCTION 


The problem of the oscillating rectangular wing in supersonic flow 
has been investigated in several papers (refs. 1 to 4) . In reference 1 
the velocity potential was derived as a power series in terms of a fre- 
quency parameter for vertical translation and pitching oscillations, and 
expressions for the aerodynamic derivatives to the third power of the 
frequency were presented. The aerodynamic derivatives obtained in ref- 
erence 2 were presented in the form of double integrals. Reference 3> 
employing the Laplace transform, presented the disturbance velocity poten- 
tial as a def ini te integral. This last result was used in reference 4 
to derive the -complete series expansion in terms of frequency for the 
disturbance velocity potential from which aerodynamic derivatives of the 
kind given in reference 1 were obtained to the seventh power of the 
frequency. 

The problem of a rectangular control surface mounted on a rectangular 
wing, which is a logical extension of the wing problem, is the subject of 
the present paper. This configuration (herein designated as a rectangular 
"wing-aileron") has been considered for steady flow in a number of papers, 
for example, reference 5* The unsteady-flow problem was considered in 
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reference 6 for a control surface hinged at its leading edge, with the 
outboard edge at the wing tip and the gap' between it and the wing at the , 

inboard edge assumed to be sealed. The results of reference 6, expressions 
for forces and moments on spanwise strips, are presented in the form of 
double integrals. A later paper, reference J, applied these results to 
obtain expressions for total aerodynamic derivatives for a wing-aileron 
combination. -- 

In the forms presented the results of references 6 and 7 are not 
particularly useful in flutter analyses employing beam theory, since 
this type of analysis calls for spanwise distributions of air forces and 
moments. In the present paper the configuration of reference 6 is con- 
sidered and, by extending the methods of references 1 and 4 to include 
a control surface, the expressions are obtained for the section lift and 
moment coefficients at any spanwise station of the wing-aileron configura- 
tion from power-series expansions of the velocity potential in terms of 
a frequency parameter. The expressions presented include terms to the 
fifth power of the frequency of oscillation. To gain some insight into 
the overall effect of aspect ratio upon the forces and moments, expres- 
sions are also given for the total forces and moments of the wing-aileron 
combination. 

The spanwise variation of the lift and moment coefficients developed d 
in the present paper are shown for a set of conditions. The section force 
and moment coefficients are applied for a single Mach number to illustrate 
the flutter characteristics of some configurations having partial-span * 

ailerons for three degrees of freedom. A comparison is provided with 
similar calculations involving two degrees of freedom. 


SYMBOLS 


A aspect ratio of wing 

b semichord of wing 

c speed of sound in undisturbed medium 


F n (x,y) 



,n-l, 


■sin" 



n 



x n- ^sin"^- 



dx 




w 


II 
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F n ,n = J ^ x n ~^sixT^ dx 


F n = F n (l,y) 


h 

h, a, 8 
la 

Jq( ) 


vertical displacement of axis of rotation, positive downward 

time derivatives of h, a, and 8, respectively 

moment of inertia of wing-aileron combination about elastic 
axis per unit span length. 

moment of inertia of aileron about x^_ (hinge axis) per unit 
span length 

Bessel function of zero order (first kind) 


k reduced frequency, bcs/V 

Li,Mi,N± components of section force and moment coefficients 
(i = 1, 2, . . . 6 as defined in eq. (26)) 

L±>Mj: components of total force and moment coefficients 

(i = .1, 2, . . . 6 as defined in eq. (33)) 


m mass of wing- aileron combination per unit span length 

M Mach number, v/c 


M a aerodynamic section moment on wing about wing axis of rotation, 

positive leading edge up 

Ma total aerodynamic moment on wing about wing axis of rotation, 

positive leading edge up 


Mg aerodynamic section moment on aileron about its hinge, positive 

leading edge up 

Mg total aerodynamic moment on aileron about its hinge, positive 

leading edge up 

2p local pressure difference between upper and lower surfaces, 

positive downward 
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aerodynamic section normal force, positive downward 
total aerodynamic force on wing, positive downward 
radius of gyration of wing-aileron combination referred to 
*o> \/lccM 2 

reduced radius of gyration of aileron referred to x-]_; v^ylp/mb 2 
static moment of aileron per unit span length referred to x-^ 
time 

flight velocity 

vertical velocity at surface of wing 

rectangular coordinates attached to wing moving in negative 
x-direction 

coordinate of wing axis of rotation 
coordinate 'of aileron hinge 


location of center of gravity of wing-aileron combination 
referred to Xq 

reduced location of center of gravity of aileron referred 
to x-^j Sp/mb 

spanwise location of inboard aileron tip 

angular displacement of wing about wing axis of rotation, 
positive leading edge up 


P = 
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6 angular displacement of aileron, measured relative to a, 

positive trailing edge down 


K = 


0 

05 


jtpb 2 /m 

density of the medium 
disturbance velocity potential 

subsidiary velocity potential associated with aileron position 


subsidiary velocity potential associated with aileron motion 


cd circular frequency of oscillation 



first natural bending frequency of wing 
first natural torsion frequency of wing 
cte natural rotational frequency of aileron 


ANALYSIS 

Outline of Problem 


Consider a thin, flat, rectangular wing-aileron configuration moving 
at a constant supersonic velocity V, as shown in figure 1. The Mach lines 
emanating from .the foremost points of the wing and aileron tips divide 
the configuration into various regions. The potentials for the various 
regions of the wing delineated by the Mach lines from the wing- leading- 
edge tips have been derived and discussed in several papers (refs. 1 to 4) . 
This paper will deal with the potentials for the various regions delineated 
by the Mach lines from the aileron leading-edge tips. A restrictive con- 
dition for the present treatment, similar to that employed for the treat- 
ment of the rectangular wing in reference 1, is that the Mach line from 
the foremost point of one aileron tip does not intersect the opposite tip. 
ahead of the trailing edge. 

A rectangular coordinate system x,y, z, moving uniformly with the 
wing in the negative x-direction, is adopted so that the xy-plane is 
coincident with the mean position of the wing, and the origin is taken 
at the intersection of the wing leading edge and the wing tip. The out- 
board edge of the aileron is coincident with the wing tip (y = 0), and 
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the inboard edge (y = yp) adjoins the wing proper. (This same configura- 
tion is treated in ref. 6.) The aileron is hinged at its leading edge, 
and the gap at its inboard edge is assumed to be sealed. 

The wing-aileron combination is considered to be performing the 
following types of small-amplithde harmonic oscillations with frequency <x> 
vertical translation h and pitching a about a spanwise axis Xq. In 
addition the aileron is allowed to perform a “rotational oso-i.l.lation about 
a spanwise axis xp. The vertical displacement at any chordwise point x 
of the wing-aileron is therefore 

z = h + a(x - x 0 ) + S(x - xp) (l) 

where the 8 term vanishes for x < xp or y > y^. The normal velocity 
at any point of the wing-aileron may then be expressed as 

w = h + Va + a(x - Xq) 4- V8 + &(x - xp) (2) 


The velocity potential for the wing-aileron configuration may be 
expressed as the sum of separate effects due to position and velocity 
of the configuration associated with the individual terms in equation (2) 
as 


0 - 0h + 0a + + 05 + 08 


( 5 ) 


The first three terms on the -right side of equation ( 3 ) apply to purely 
wing oscillations and were treated in power-series form in reference 1 
to the third power of the frequency and in reference 4 to the seventh 
power. The last two terms on the right side apply to control-surface 
oscillations and will be treated in the following section. 


Aileron Velocity Potentials 

Consider the flow field of the aileron, for which and ^ 

are the velocity potentials. The field is divided into three regions 
(fig. l) : region I behind the outboard Mach line, region III within 

the inboard Mach cone, and region II which is the remaining part of the 
control surface. Expressions for the velocity potential in these various 
regions will now be obtained. 
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Regions I and U ._ In accordance with linearized theory the velocity 
potential for regions I and II can be obtained directly from the velocity 
potential of a semi-infinite rectangular wing. This is because the aileron 
is not affected by the presence of the wing when the wing is at constant 
zero angle of attack and is therefore essentially a wing itself. • The 
boundary- value problem for the semi-infinite rectangular wing was solved 
by means of the Laplace transform in reference 3 < The transform of the 
velocity potential was obtained in the form 


where 



w [l _ 2_ f* /2 e -|^y sec 2 e 

pU " itn J 0 



(4) 


Vi = 



a£ 

C 2 


P 



and where the bar above w, 0, and x denotes the Laplace transformation 
of these terms with respect to x. 

Employing the 5 and 5 components of equation (2) in equation (4) 
and applying the inverse Laplace transform (pairs 55 and 84, pp. 298-299, 
ref. 8) result in the following expressions for the aileron potentials $5 

and in regions I and II: 



e- i5K dx 
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Note that the upper limit of the second Integral, -which is within 
the brackets of equation ( 5 ), will vanish when x equals (3y. This 
indicates that the second integral represents the correction to the 
"purely supersonic" disturbance potential needed to compensate for the 
flow around the wing or aileron tip. The first term is the purely super- 
sonic disturbance potential obtained in reference 9 * 

The indicated integrations in equations ( 5 ) and. (6) do not appear 
to be obtainable in terms of known functions. For purposes of evaluation, 
however, these equations can be expanded in a power series and terms can 
be retained to any desired power of the frequency parameter cb. For 
region II, the region of "purely supersonic' flow wherein x ;> 3y, equa- 
tions ( 5 ) and (6), expanded to aP, yield 


05 = * -|r f l< x 2) (7) 


V 
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0s = - | f 2(x 2 ) 


( 8 ) 


■where 


fl(x) 


x 


lose 2 

2 


+ 3) 4- S*sW + 5) 
lasr hSusr 


+ ^r-(8p 4 + 4op 2 + 55) 

960M 4 


- 9^ E ^ _(8(3^ + 56(3 2 + 63) (9) 

576014 ^ 


f 2 (x) = 



( 1 °) 


and x 2 = x - x^. For region I, the region of "mixed supersonic" flow, 
equations (5) and ( 6 ) yield 


$5 


2b Y5 

0 


Al(x 2 ,y) 


( 11 ) 


$5 = - — p M x 2>y) 


( 12 ) 
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where 


Ai(x,y) = 


155F 2 - ^-(xF 2 + 5x^2 + (2P 2 “ 1)J 

ar ' 7 iarL J 


of 


48m 


r 4L 


x?(4p 2 + 5)F 2 - 3xFi| + 2$%?^ 


lap 


960M 


4L 


5x^(4p 2 + 7)Fr 


30x% 4 + (8p 4 - 4p 2 + 5)Fg 


A 2 (x,y) = f \ 1 (x,y)dx 
J n 


(13) 


(14) 


The functions fj^x), f 2 (x), A-j^x^y), and A 2 (x,y) were given to the 

seventh power of 5> in reference 7 and are repeated here to oP. In 
equations (15) and (l4) 


and 


F n ( x >y) = 



x n-1 sin -1 ^^ dx 



m 


^-(x m+1 - Fnrt-n+l) 


Further information regarding these F n functions can be found in ref- 
erences 1 and 4. 

The above equations could have been obtained from the equations of 
references 1 and 4 by substituting 5 for a, 0 for Xq, and x 2 

for x. In other words, to obtain $5 and the aileron in super- 

sonic flow may be considered as a rectangular wing rotating about its 
leading edge insofar as regions I and II are concerned. 
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Region III .- The problem of determining the -velocity potential for 
region III may be complicated by the existence of flow about the inboard 
tip of the aileron. The present treatment of the problem avoids this 
complication by assuming the’ gap between aileron and wing to be sealed 
as was done in reference 6. 

The velocity potential for region III can be obtained from refer- 
ence 9* In reference 9 the potential of a "mixed supersonic" region at 
the tip of a rectangular wing in supersonic flow was considered. If 
this region is assumed to be part of a "purely supersonic" region, an 
approximation for the velocity potential in region III can be obtained. 
Although the expression so obtained was an approximation for the velocity 
potential in the tip region then under consideration, it is precisely the 
potential for region III of the current problem. This expression, as 
applied to the present case, is 


4 = — f 2 (VS + £6)e -i “( x 2-£) f cos 

fa J 0 Jo 


g(x 2 - £)sin 8 


d6 d| 


1 r x 2-PY2 

Pit J 0 


(Vb + ib)e-^2^) f 


cos 


-1 &Z 


X2-I 


cos 


— (xp - I) Bin 8 d0 d£ 
M c 


•]' 


- / 2 (V6 + !6)e“ la) ( x 2“^ J 

P4o ° 




1_ r 2 Py2 ^ + ^ e -io)(x 2 -£) 

faJo 


. cos 


d| - 

-1 ^2 

X 2~ ^ cos 




|)sin 0 


de d£ 


( 15 ) 


where 

72. = vi - y 


(This expression for the potential applies to those points of. region III 
that are on the wing as well as those on the control surface.) 
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The first term on the right-hand side of equation (15) is the sum 
of the potentials $6 and $5 of the "purely supersonic" region, 

region II, given in expanded form in equations (7) and (8). The calcula- 
tion, therefore, of only the second term is necessary. After expansion 
of equation (15) in a power series to the fifth power of 5>, the velocity 
potentials for region III can be expressed as 


V6 


+ - sin -1 — jfqtxg) + — cosh -1 f 5 (y 2 ) - 

\ 2 * x 2 / « P|y 2 | 




^\l*2 2 - ^2 ^(x 2 ,y 2 ) 


(b + ~ sin -1 — )f 2 (x 2 ) + — cosh" 1 — — 

\ 2 * * p|y 2 | 


(16) 


^5^2) + %(y2) 


^~\j x 2 2 “ p2y 2 2 A ^( x 2»y2) [ 


(17) 


where 


fj(y2) = 


%(y2) 


pV 2 2 ^ P 8 y 2 ^ 

— _ — + — t — 

2 4 

12M 320M 



8 om £ 


(18) 


6 


( 19 ) 
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and 

^(x 2 ,Y 2 ) 


\(X 2 ^ 2 ) 


ch 2 *. p icb^Xp 2 io^p^p 2 

^ + ?(P 2 + 3) - - : — ~(2p 2 + IT) ^“(2P 2 “ 1) 


2 ’ 2 

lar 


itoM c 


7 ar 


— — ( 2 p^ + top 2 + 45) - ^=- (3P k - 10 ) + 

9601 # 96 om^ 


ia?Xp il ' j. P ia^Xp^^p 2 2 2 

2 / QoA , onf^o^ j- )■ N j. 2 =— (l6S^ - 8(3 - 219) + 


28800M^ 


-(8p + 296p^ + 453) + 


13200M )+ 


- 4*2 + S 

10800M 4 


( 20 ) 


-22 -222 

iciKp OJ Xo _ 05 p y P 

1 + -jS + — A(2p2 + 5) - -f_|-(2p2 + 1) - 

2 5 ■ 4 a * 2 2 tor 


iai?x 2 ^ 

720M 2 


(to 2 + 25) - 


io^x^ 2 y 2 2 


(lip 2 - 10) - 


720M C 


-4 4 -4. 2 n 2„ 2 

05 Xr 


57601 ^ 


^to 1 * - + 5to 2 + 61) - 2 P — (np^ - 8p 2 + 43) + 


86toM H 


°° P y 2 -(8p^ + 4p 2 - 1) 
2l60M^ 


( 21 ) 
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Other regions .- Two other types of regions are possible. One 
region is formed "when the Mach line from the inboard edge of one aileron 
intersects the inboard edge of the aileron on the opposite side of the 
wing. This type is not of practical interest, however, since the inter- 
ference effects of the fuselage between the ailerons would be of far 
greater importance. The other type of region, formed when the Mach lineB 
from the tips of one aileron intersect on the aileron, presents no diffi- 
culties. The potential in this type region is simply the sum of the 
expressions for the potential in regions I and III minus the potential 
for region II. This parallels directly the treatment of region Tj in 

reference 1. 

Expressions for the velocity potentials due to aileron motion have 
been provided for all regions in the field, of flow of ' the aileron ( see 
eqs . (7), (8), (ll), (12), (l6), and (17))* The aforementioned equations 
are in terms of known integrable functions and, as subsequently shown 
in the following section, yield force and. moment coefficients for a 
chordwise strip that are valid for sufficiently small values of the 
parameter cd. 


Section Forces and Moments 

Expressions for section force and moment coefficients for any 
station along the span are derived. In deriving these expressions the 
variables x, y, Xq, and xq are employed as nondimens ional quantities 

obtained by dividing the variables previously used by the wing chord 2b. 
The local, perturbation pressure difference between the upper and lower 
surfaces of the wing may be written as 


-2p(^ + X^ 

\8t 2b 8x > 


( 22 ) 


where p is the density of the medium. The section force (positive 
downward) on the wing is therefore 


P = 2b / Zp dx 
J 0 


( 23 ) 
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The section moments (positive leading edge up) acting on the wing section 
about the wing axis of rotation x Q and on the aileron section about the 

hinge point x^_ are, respectively. 


Ma = 



Xq)Zs) dx 


(24) 


Me 


= 4b 2 (x - x-j_)/3? dx 


(25) 


After substituting equations (3) and (22) into equations (23), (24), 
and (25) and performing the indicated integrations, the results may be 
written as 


P = -4pbV 2 k 2 


— (Lq + 1X2) + ot(x^ + 1x4) + s(x^ + 1x5) 


= -pV^bc, 


M, 


a 


= -4pb 2 V 2 k 2 


^(M-l + iM 2 ) + a(M 3 + 31%) + 8(145 + 3Mg) 


= 2pV2b2c m 


(2 6) 


Mg = -4pb 2 V 2 k 2 


^(N-l + iN 2 ) + a(Nj + il%) + S(N 5 + iNg) 


= 2(1 - x-^pV^Si 


The terms Xj_, M-j_, and Nj_ (where i = 1, 2, . . .6) are the in-phase 

and out-of-phase components of the force and moment coefficients; the 
odd-number subscripts are the in-phase components and the even-number 
subscripts are the out-of-phase components. Terms with subscripts 1 
and 2 are associated only with vertical translations of the wing. 

Terms with subscripts 3 and 4 are associated with angular position 
and rotation of the wing about x = Xq. Indices 5 and. 6 are associ- 
ated with the angular position and rotation of the control surface 
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about x = -Xq. As may be seen from equation (26) , the forms L^, Mj_, 

and Nq, which have been used in studies of unsteady flow, can be readily 
converted to the conventional section coefficients for lift c i, pitching 

moment c^, and aileron hinge moment c^. 

Expressions for the coefficients of equation (26) differ for the 
various regions of the wing and aileron. Coefficients associated with 
wing motion, namely Lq, L 2 , L 3 , L 4 , Mi, M 2 , - Mj, and M 4 , were 

given to the third power of the frequency parameter tb in reference 1 
and were extended to the seventh power in reference 4. For completeness, 
however, these expressions are repeated to the fifth power of in the 
appendix to this paper. 


The remaining coefficients can be divided into two groups: (l) those 

which show the effect of wing motions h and a upon the aileron 
moment M 5 , namely Nq, N 2 , N^, and 3%, and (2) the aileron section 

coefficients due to the aileron potentials and namely L 5 , Lg, 

M^, Mg, N 5 , and Ng. Expressions for the coefficients in each of these 

groups for the various regions of the wing-aileron combination will now 
be given. 


Aileron coefficients Nq, N2, Nj, due to wing motion.- For 

the section coefficients Nq, N2, N3, and N^, the aileron is divided 

into three regions (fig. 2). Region 1 includes all sections outboard 
of the point where the wing-tip Mach line crosses the aileron hinge, 
region 3 comprises all sections inboard of the point where this Mach 
line crosses the trailing edge, and region 2 includes all intermediate 
sections. The coefficients can be expressed as follows:' 


N i = 


N n 


*1 


- 2 x-iL 


1^1 


N 2 = N2* + M 2 1 - 2xqL2* 

N 5 = N 5 ' - 2 x 0 Nq 


( 27 ) 


N4 = - 2xoN 2 I 


where the coefficients Mq 1 and Lq' (i = 1, 2, 3> 4) 
wing coefficients given in the appendix for Xq = 0 . 




refer to the 
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For region 1 where 0 < y < X;L /p, the coefficients of equation (27) 

SntS in for °°f flclmts ^7 

and are given as the foilL^ e^sl^! *" tte rSgl ° n 0 < * < Vp 


N l' =_ i^ + 1)F 2 " + ^- + 2 %" + 


0 * 

aA 2 


m 2 ’ = ± 
* 0 * 


'[ X l^ 2 " “ Sxi 2 ^” + (8p 4 + - D^" - 0 2 (5p2 + 4)F 5 'jJ 

^ x i F i' ~ f 2^ + ^[ (2p2 " - 6 v\? 3 " + (4p 2 + + 

--^-[5x1 V 2 " - 2x^(6^^ - 4p 2 + 3^" + top^F^' - (sk0^ + 8p 2 _ x) F g”J 


N 3 ' = M 3 ' . 2x^3' + it 

“"Pit 


i( Xl Fi" - F 2 ") - ^ Fx" + 


x 2 If ^ 

"i: + 3P 2 - 1) - Xl (6p2 + 5) + ~(20p^ + 2ip 2 + 3) + 


P £ 


30" 


o of 

^[*lV(0 2 + 5) + + £0^4 _ 5p2 + 3) + 

^^.■^"(^P 2 + 13) - (5 6p^ + 64p^ + up 2 _ i)Fg" 

%' = M4 1 - 2x31.4 ’ + jjj- |^KV - ^ F 2 M (30 2 + 1) + (2P 2 + 1)F3 -j + 

~[^2 n (20 4 - P 2 +1) - lax^sy + * A -(aop^ + 7p 2 . X) _ 

2 p 2 F 5 "( 5 p 2 + 3 )J _ ^~[3bci5F 2 "(p 2 + 7 ) + 10 Xl 4V r (8p 6 - 4 p 4 + 

9op 

30 2 - 3) - 36op^xx 2 F3" + 3xxFg"(l68p^ + 64p lf - up 2 + 3) _ 

^ k Fj”(28fi 2 + 19)] 


(28) 
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where 


v- r 1 * 1 " 1 
J 0 

For region 2 where Xq/p < y < l/p, the terms in equation (27) are 

calculated by again using the expressions for the wing coefficients pre- 
sented in the appendix (eqs. (Al) and (A2)) for the region 0 < y < l/p 
and yield the following expressions : 


■1,/W 


dx 


N-i‘ = 


_ 4p 2 + 3 
3P 5 15 p^ 






x -, 2 M 2 kx n ^ M^k^x-i ^ 

' - - + -(4f)2 + 7) 

0 k 2p5 36p9 


Xq 2 X l^ o M^Xq 6 ), „ 

N 3 ' = M 5 ' - 2 xqL 3 ‘ + — — (p 2 + 3) + — + 35P 2 + 35) 

pk 2 6p 5 180 P 9 


) (29) 




1*4* =.!%' - 2xqlV + — — (p 2 - 1) + — (p 2 + 5) 

15P 7 


i_C4p^ + 49P 2 + 63) 

630P 11 


For region 3 where l/p < y, the terms of equation (27) are obtained 
by using the wing coefficients in the appendix (eqs. (A3) and (A^t-) ) for 
the region y > l/p and equation (29) « 
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Aileron coefficients Lg, Lg, Mg, Mg, Ng, Ng due to aileron 

motion . - The expressions for the aileron coefficients Lg, Lg, Mg, 
Mg, Ng, and Ng depend upon which region (see fig. l) contains the 
trailing edge of the section. 

For a section partly in region I, that is 0 < y < Xg/p (where 
xg = 1 - Xj_) , the coefficients are 


- — 

“ p* 


- gV' + + » 2 - « - ^-^e 2 + 5) + + 5) + - 

XjV'Ctop 6 +■ aop 4 - 50 2 + 3) + 2p 4 7 5 "'<15P 2 + 13)] 


*-fc 


ac^ 1 ' • _ * 2 e -(jpg + i)l + ' , (2s 1 *' - p 2 + i) - sftp'Sgiy ' ' + " (a>p 4 + 7P a - i)j - 

P J 3p 

^ 6 [^ 3 1| T g "'(P a + 7) + 10*3%' "(Ep 6 - 4p 4 + as 2 - 3) - Stop^Fj"' + J 6 ’"(l68p 6 + 64p 4 - lip 2 + 3)]| 



XI III 
*2 

p* 

l* 2 




I 2 "' - * tX ? 3Fl ' + X ? (6p 4 + 3P 2 - 1) - ^t-^aop 4 + Zip 2 + 3) + ^(kp 2 + 5)^2 ‘ 1 ' + 

| k 2 3 p 4 ^ 12P 8L 

3S» 6 x^'" - 2x3 S F'i t " , (4op 6 + asp 4 - 5P 2 + 3) + ?6'"(56p 6 + ftp 1 *' + UP 2 - 1)] j + 2 (*i - *0)^5 

- (sp 2 + l)!."'] + 2S£|’&Jp 2 , ' , (ap*- - p 2 + 1) - 12pWV'' + 2P% , 5 ”'(3P 2 + 3)} - 


. (30) 


2^[3x 3 5 F a ,,, (p a + 7) + 5 x 3 ^'" (Sp 6 - V* + 3P £ - 5) - sop^ij^j'" + 2 p 4 T 7 ,,, ( 28 p 2 + 19)] 
45P 10 


+ 2(lj_ - 


lg “ M3 “ 2(x]_ - XqJIj 


»6 - «6 - 2(*L - ^ 0)^6 
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where 


, x 


P. 


n 


^ x n “^- sin dx 
x 


X? 


For a section entirely within region II, that is — < y < y x 


these coefficients become 

x-t x-z? k%^x z 3 

L 5 = _^p2 + 3) + — 

pk 3P^ 6op' 


s' 

2-(4p^ + 35P 2 + 35) 


^6 






Xj^ 0 km,“3s , 0 „ K-'jyi ‘x*“ , 0 

= — (r 2 _ 1) + 3 -(p 2 + 5 )• ~( 4 p 1, ‘ + 49 p^ + 63) 

■* n - JJ. 


kp- 


6p? 


l8op 


x 2 x 4 j c ^ I 2 x 6 

Me; = -JL - ~(p 2 + 3) + y-j~-(4p^" + 55P 2 + 35) + 2(x X - Xq)L 5 


5 pk 2 2p 5 


36p9 


^ = ! i£( P 2 . 1 ) + !?^£( P 2 + 5 ) - ^ + 63 ) + 


3kp' 


15P 1 


105P 


11 


2(*i - x 0 )Lg 


Npj ss - 2(x x - Xo)Lij 



f (3D 


N6 = M6 - 2(x x - Xq)L6 
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For a section partly within region III, that is y, - < y < y-i + -2, 

1 0 1 P 

these coefficients become 


L R = 


x 3 x 3 3 


--T o 

— — r(P 2 + 3) + — -(4p^ + 35p 2 + 35 ) 

pk 2 3p5 6op9 


1 + 1 sin" 1 + 

2 it Xz / 


y 2 x 3 

— cosh - - 1 - 

3C Py 2: 


1 - 2Xj 2 . ^_(Sp2 , + 5, + 

k d 3 


ai2k 2 y 2 2 x ^2 


M 2 k 2 y 2 ^‘ 


20 


(15P 2 + 13) 


♦ XJ r\ft - * v 


—(12^4 + 4p2 - 3) + 

3^ 


M^k^x, 2 x- j, 0 l^k 2 ^ 2 ^ , 

-2— (lop b - 22p 4 + 55p 2 + 135) L§-( 265 P 6 + 83p 4 - 50p 2 + 30) 


l8op 


l8op c 


( 32 a) 


l 6 = 


kM^x-z^ 


k^M^Xz^ 


— -(P 2 - 1) + ~ ■ (P 2 + 5) frC^P 4 + ^9P 2 + 63) (i + 

(kp2 6 r 1 T Ann H 1 2 


6pT 


l8op 


I sin-1 522) + 22S oosh -l 2, 

It Xt It R- 


Py 2 


4 Si 2 ky 2 2 ^Sl^y^ 


r i\F? - ■> V 


5P 




L“ ' 15^ ' up2 + 17) - 


liM2y 2 2 / k o M^k^x, 1 *- , 

~-(20p 4 + TP 2 - 1) + 2_(8p6 _ tfpk + 159p 2 + 453) + 


9P 

li^^x^Syg 2 

l8op 8 


900P 


10 


< 126 p 8 - 17 2p^ + 203P 2 - 219) + 


34^k5y 2 4 

675P 6 


(l68p 8 + 64£^ - up 2 + 3) 


( 32 b) 
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Mr = 


*£ - + 3 ) + + 350 2 + M 


pk 2 2p5 


eoeb- -l_|l 

3k Py2 


,-l x 3 


36p 9 


M 2 k 2 y 2 2 ) t Y2 


1 + 1 ein-l !£\ - 

2 k 


TvT3 


far 


\ + ^l- ^ (88p i, • + 


_k 2 l8p^ 


33p 2 - 21) - ^-(20p it - + 2ip 2 + 
9P 2 


3 ), + M k x | _(i6p g - 36^ + 8ip 2 + 209) - 
l8op 8 


M2 k X 3 2y 2 (),A«p6 + i72p4 - 97 P 2 + 47) + “=-{560 8 + 6^0 ^ 


270 p c 


lip 2 - 1) 


1350 


+ 2(x]_ - Xq)L^ 


(32c) 


M6 = 


4x^3 0 4M2kx,5 M 1 +k3x 3 T 

* n d-(p2 + 5) 2-(4p4 + 49p 2 + 63) 


^(0 2 - l) + 

3k0^ 150 ' 

i sin -1 — ^ + — cosh -1 

k x, i k 0y 2 


1050 


11 


1 + 

2 


3 


W 2 + 1) - + 


3k 


A4^ (2ep2 + 19) 


126 


£2^* j 2 - fV 


i_(2p 2 + 1) - 


30^ 


- 76 2 + 10) - 1^4562 + 5) + - ^P 4 + 

9450 1° 


150 


50 £ 


2310 2 + 77T) + 


lj^3^ 2 -(i33p 6 - 1850^ + 2100 2 - 210) + 

n)i Kn 8 


9450 


M^k3y p4 

^-(28p 2 + 19) 

1260 2 


(32d) 


+ 2(x x - x 0 )Lg 
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1*6 = Mg 


2(xi - Xo)L 5 

( 32 e) 

M° 

1 

CVi 

(32f) 


Equations (32) apply to all sections of region III - those sections on 
the wing as well as those on the aileron. 

The case of a section in a region within the Mach cones from, both 

x 3 x 3 

aileron- leading-edge tips y-, - — < y < — , which is formed when the 

x P p 

Mach lines intersect on the aileron, can be treated with the expressions 
provided in this section. The coefficients for such a section are obtained 
by adding the corresponding expressions of equations (30) and ( 32 ) and 
subtracting the pertinent expression of equations ( 31 ) • 


Total Forces and Moments 

In order to gain some insight into overall effects of aileron loca- 
tion and aspect ratio, expressions for total forces and moments are pro- 
vided. The total lift, total wing moment, and total aileron moment, 
respectively, can be written as 


P = -8pb 2 V 2 k 2 A 


^(Lq + fL 2 ) + a(l>5 + iL4) + ^ 5(I 5 + iL6) 


M a = -8pb^V 2 k 2 A 


-(Mq + 1M 2 ) + a(% + ii%) + 5(M 5 + iMg) 

b bA 


yi 


> (33) 


ii s = -8pbVk 2 & 

b 


— (Nq + iif 2 ) + a(% + ii%) + 5(% + iNg) 
b 


where A is the wing aspect ratio. Expressions for the coefficients Lq, 
L 2 , Lj, l k , Mq, M 2 , Mj, and M^ in equations (33) are given to the 

seventh power of the frequency in reference For completeness they are 
provided to the fifth power of the frequency in the appendix. 
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The other coefficients associated with wing motion are 


2 x - j 3 

M 2 k 2 x 1 3 

pi a 2 4 . 

b 

k 

3p3 

15P 7 

37 i 

2p4 


x > _ M 2 ^! 6 . 

Ji _( p 2 + 2 ) + 

i 5 po 


8p 2 + 8) 


+ M ^' - 2 ( x ! - Xq )^' 


( 34 a ) 


Xn 2 M ^ kx -,4 M 4 k 3 x 1 ^ 

+ T^(4p 2 + 7) 


pk 


2p- 


56p 9 


b 


x-t 3 M^kx-i 5 
-i ^-( p 2 + 4 ) + 

p 2 k 5 P 6 


2 M 4 k 3 x ]_7 


105P 


10 


(p4 + I 2 p 2 + l 6 ) 


+ M2 1 - 2 ( x ! - Xq )^' 


(34b) 


X n 2 Xn 4 _ M 2 k 2 x - 1 ^ / i , o 

_i L_(p2 + 3) + — (4p4 + 35P 2 + 35) 

pk 2 6p5 I80p9 


b 

^1 


x r 


Xn 5 ■ 

-~(3P 2 + 4) + 


+ ao P 2 + i «) 


ao 


p 2 k 2 

+ Mj' - 2(x x - x D ) (£5' + 


5P° 105P- 

2xoLi') - ^ (N X - Mi') 


^(pa . 1 ) + Mgl ?-£(pg + 5 ) - + 49P 2 + 65) 

3 v i ~7 6 J 0 P 11 


(34c) 


jp-'k 15 P 

x / J + kM ^ 6 


b 

3 yi 


m^ 5 x 8 

(p 2 + 2 ) + ri-(3P^ + l 6 p 2 + 16 ) 


»12 


p4k 15 P 8 105 P 

M4' - 2(x x - Xq)(iV + 2 xoL 2 ') - 2 xq(N 2 - M 2 ') 


(34d) 
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The primed, coefficients employed in equations (34) are obtained by 
replacing the wing aspect ratio A vith y x /b in the corresponding 

unpri me d coefficients in equations (A5) and (a 6) of the appendix. The 
remaining coefficients, which are associated with the total forces and 
moments on the aileron, are 


L 5 = 


x, k %^x^5 

J L(p2 + 3) + + 353 2 + 55) 

pk 2 33 5 6op? 


+ 4) + p* + 20P 2 + x6) 


b / 2' 

— (1 + — 

v 


p 2 k 2 


53 c 


453 10 


L 6 = 


^£(p2 _ 1) + m2x 5 . , (p2 + 5 ) _ ^ -(4p 7 + 493 2 + 63) 

?k 6p 7 ^ 


-S-fi + s) 

M *y 




p^k 

2- _ ^(p 2 + 3) + 
pk 2 2p5 


acj3 _ to^£ (|3 2 + 2 ) + + X6p2 + X6) 


2 4 

Xz^ X ^ 


53 8 

k% 2 x* 6 


56p 9 


103P 12 

2-(4p^ + 353 2 + 55) 


b (l + 2 ^ 


4 x* 5 4k% 2 x 5 7 p 

— 1-(3P 2 + 4) + — (53 4 + 20p 2 + 16) 


53 


553 


10 


5y x \ 

+ 2(x x - Xo)l 5 




p 2 k 2 


M6 = 


- X) + ?^4 2 ♦ 5) - + & * 63) 


5kp3 


153 7 


1053 


ll 


b L 2> 

— 1 + - 
5yn \ 


3x^4 4kM 2 xx 6 M^k5x 3 8 . 2 

5 — (p 2 + 2) + ^-(53^ + 3-6p 2 + 16) 


kp i 


53 


8 


153 


12 


2(x x - Xo)L6 


- 2(x x - Xq)L 5 
K(J = M6 - 2 ( x l " *0)^6 
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In equations (35) the expression -within' the first set of brackets in 
each equation represents the component of the coefficient due to two- 
dimensional flow, and the second term provides the corrections for edge 
effects. The correction for the inboard or sealed edge is 2/jt times 
that for the outboard edge of the aileron. This serves as a check. on 
the calculations, inasmuch as it is in accord with the observations of 
references' 6 and 9* Equations (35) do not include the contribution of 
the part of region III that is off the aileron. If this contribution 
were included, thereby providing overall total coefficients for the wing- 
aileron combination; the correction for the inboard edge would vanish. 

For an aileron whose outboard as well as inboard edge- can be considered 


sealed, the expression 1 + — in equations (35) should be replaced 

jt 

by 4/rt. 


CALCULATIONS AND DISCUSSION 
Spanwise Aerodynamic Effects 


To give some indication as to the general nature of the spanwise 
distribution of the different components of the lift and moment coeffi- 
cients introduced by the presence of a control surface, these components 
were evaluated for a given set of conditions at different spanwise loca- 
tions. This spanwise variation is illustrated in. figure 3 for an aileron, 
with a chord one-half of the wing chord and a span one-fourth of the wing 
span, mounted adjacent to the wing tip on a wing of aspect ratio 4 at a 
Mach number of 1.3 and a reduced frequency k of l/9- Figure 3(a) pre- 
sents the spanwise distribution of the components which contribute to 
the wing lift force and moment, and figure 3(b) presents the distribution 
of the components contributing to the aileron hinge moment. Two- 
dimensional values of the coefficients are represented by dashed lines. 

For the aforementioned set of conditions, the Mach lines from the aileron 
tips intersect on the aileron. (For stations partly within the region 
formed by these intersecting Mach lines, see the discussion following 
eqs. ( 32 )). 


Application to Flutter 

The application -of the section force and moment coefficients devel- 
oped in this paper to some particular examples of the flutter of a 
finite-span wing with partial-span control surfaces is considered of 
interest. A broad systematic study, however, such as that of reference 10 
which employed two-dimensional aerodynamic theory, is beyond the scope of 
the present paper. 
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The particular examples deal with a rectangular wing-aileron of 
aspect ratio 4 at a Mach number of 1.3. This wing, without an aileron, 
has been analyzed for flutter and the flutter characteristics have been 
calculated and compared with experiment at M = 1.3 in reference 11, 
where two-dimensional flow coefficients were employed, and in refer- 
ence 12, where finite-span coefficients were employed (refer to model A-l 
results of the reference papers). The following pertinent parameters of 
reference 11 are used: 


Aspect ratio, A 4.0 

Elastic-axis location, Xq, referred to 2b 0.413 

Center-of-gravity location of wing-aileron combination^ x&, 

measured from elastic axis in semichords O.I56 

Radius of gyration of wing-aileron combination, r a = y^/mb 2 , 

referred to Xq 0.28 

First natural frequency in bending, <%, radians/sec 838 

First natural frequency in torsion, o^, radians /sec 1,747 

Density parameter, i/k = m/jtpb^ 64-9 


The present study considers the wing described above to have an 
aileron, hinged about its leading edge, with a chord one -half the wing 
chord and a span one-fourth the wing span. The aileron is considered 
rigid but spring restrained at the hinge. The effect upon flutter of 
varying the spanwise location of this aileron is examined for two mass- 
balance conditions. The aileron parameters assumed for this study are 


yi/ b 2 

Aileron-hinge location, x^, referred to 2b 0.5 

Reduced aileron center-of-gravity location, Xp, measured from 

hinge in wing semichords 0.2 and 0 

Reduced radius of gyration of aileron, rp = |/lp/mb 2 , 

referred to x^_ 0.125 and 0.1 


Flutter analyses for three degrees of freedom in three-dimensional 
flow . - Three degrees of freedom - namely, wing bending, wing torsion, 
and aileron rotation - are involved in the present study. By the method 
given in chapter IX of reference 13, the flutter determinantal equation 
for the three degrees of freedom involved may be written as 
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A^ Acta 
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■where the determinant elements are 

A*h - [i - fe) 2 ] - s * fa + 

Ata - h Z a ay - \ *f 0 ( L 3 + ay 

Ahp - Xp/^z h% ay - t «/^ 2 (l5 + “6>Vp ^ 

A ob = Xa,/ Zo^h ay - £ * (% + 1M 2> Z a Z h 

V* - r a *[i - (h) 2 ] / Vay - i «/ 0 V + 

V . jy + 2xp( xi - X.)] f\ft ay - £ ■ /^ 2 (“5 + ay 

Ap h = xp f\ z h ay - i k/JW + 4y 

V = [*/ ♦ **(*! - ^ ffa a ay - | + “JV. * 

- *4 - ftflfa** - * + 47 




NA.CA tst 3 644 


29 


and where Zq and Zp represent, respectively, the inboard and outboard 
edges of the aileron, Z^, Z<x, Zp represent, respectively, the mode 

shapes in wing bending, wing torsion, and aileron rotation, and the coeffi- 
cients Lq, Mp, Nq are defined in equation (2 6 ) . In this analysis the 

wing modes Zi^ and Zq, were assumed to be the uncoupled first bending and 
first torsion modes of a uniform cantilever beam, while the aileron mode 
was assumed to be unity. The aerodynamic integrals required in these equa- 
tions are obtained by numerical integration employing the basic method shown 
in appendix B of reference 4. 

After substitution of these integrals and parameters into the deter- 
minant elements, the equations are solved, as in chapter XIII of refer- 
ence 13 , at several values of the reduced frequency k for the . flutter 
frequency to and aileron frequency cqg. The results are presented as 

plots of the parameter b<%/e against the natural frequency ratio 
in figures 4- and 5 . The parameter bc%/c is employed here as in refer- 
ence 10, since It is preferable in a compressibility study rather than 
the velocity parameter normally used in incompressible flow. The choice 
of <% rather than c% in this parameter is purely arbitrary since in 

this analysis <%/(% is constant. 

Figure 4- illustrates the results of the flutter analysis for the 
mass-balance condition defined by the parameters Xp = 0.2 and r ^2 = 0.125 

implying a relatively rearward aileron center-of -gravity location. The 
analysis was performed with the aileron in three spanwise locations on 
the wing, as shown in figure 4. In figure 4(a) the aileron is in the 
Inboard position adjacent to the wing center line. In figure 4(b) the 
aileron is located midway between the wing center line and the wing tip 
and In figure 4(c), adjacent to the wing tip. 

Two curves appear on these plots and regions of instability relative 
to each curve are shown. The lower curve extends through the entire 
range of the frequency ratio ug/o&, and is associated with a mode that 

is predominantly a combination of wing bending and wing torsion.. As ug 

becomes Infinite, the case of coupled wing bending torsion Is approached; 
the value of boag/c approached asymptotically is Indicated by the short- 
dashed line and corresponds to the result from reference 12. As erg 
decreases the value of bo^/c increases slightly. At a sufficiently low 
value of a point is reached where the region of instability associated 

with the upper curve appears. Since the stiffness bo^/ c associated with 

this upper curve is higher than that associated with the lower curve, the 
mode associated with the upper curve is the critical one for low values 
of the frequency parameter ug/c%. 
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In figure 4 parts of the curves are estimated as Indicated by the 
long-dashed lines. -Over the estimated portion of the curve the reduced- 
frequency parameter k is greater than 0.25 and beyond this value the 
aerodynamic derivatives developed in the present paper become inadequate. 

To provide some indication of the effects of mass balance, the con- 
figurations of figures 4(b) and 4(c) are again analyzed for the mass- 
balance condition defined by the parameters _x^ = 0 and r^2 = 0.100, 

which places the aileron center of gravity at the hinge. These results 
are shown in figure 5 where only the cases of the aileron at the midspan 
location (fig. 5( a )) and at the outboard location (fig. 5(b)) are 
considered. 

For figure 5(a) the maximum value of 02 g/c% for the upper curve 

is unknown and hence this curve is not closed. A comparison of fig- 
ures 5(b) (balanced aileron) and 4(c) (unbalanced aileron) shows that 
the upper curve is extended to a much higher value of ■ in fig- 

ure 5(b), and indicates that a stiffer (higher frequency) aileron is 
required for stability. This finding iB in accordance with an observa- 
tion of r.eference--10 that mass balance may have an adverse effect for 
low supersonic Mach numbers. 

Two -degree-of -freedom subcases .- It is often deemed practical in 
examining the problem of flutter for a multidegree-of- freedom configura- 
tion to consider pertinent subcases under the assumption that the char- 
acteristics of the solution to the multidegfee-of-freedom problem will 
manifest themselves in the solutions to the pertinent subcases at some 
saving in labor. To test this assumption the subcases of the present 
three -degree -of -freedom problem which involve aileron motion - that Is, 
the wing-bending — aileron-rotation case and the wing-torsion — aileron- 
rotation case - have been considered. 

The flutter determinants are obtained for the bending- aileron case 
by eliminating the row and column associated with wing torsion from the 
three -degree-of -freedom flutter determinants employed in the previous 
sections. Similarly for the tors ion- aileron case, the row and column 
associated with wing bending are eliminated. The results of this study 
are shown in figures 6 and J. 

The results obtained by using the finite-span coefficients of the 
present paper for both mass-balance conditions are presented for the 
wing-bending — aileron subcase in figure 6 and for the wing-torsion — 
aileron subcase in figure J. The curves in figure 6 have the same form 
as the upper curves in figures 4 and 5 and, about the same critical 
ranges of bcj^/c. For the aileron at an outboard location a comparison 
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of figure 6(c) with figures 4(c) and 5(b) provides a somewhat interesting 
observation; namely, that, when only two degrees of freedom are consid- 
ered, the critical frequency ratio for the first mass-balance condition 
is increased, whereas for the second or mass-balanced condition the 
critical frequency seems to be slightly decreased. 

In figure 7 are presented the results for the tors ion- aileron case. 

No region of instability was detected for the aileron in the inboard case. 
Therefore, results are presented only for the midspan location (fig. 7(a)) 
and outboard location (fig. 7(b)). The curves in these figures appear 
to correspond in form to the lower curves in figures 4 and 5* Whereas 
mass balancing the aileron did not appear to materially affect this 
stability boundary when three degrees of freedom were considered, for 
the mode under consideration a mass-balance control surface would seem 
detrimental from a flutter standpoint when only two degrees of freedom 
are considered. 


CONCLUDING REMARKS 


The present paper has extended the treatment of NACA Rep. 1028 
and NACA TN 3076 for rectangular wings in supersonic flow to include a 
third degree of freedom, namely partial-span aileron rotation. Expres- 
sions for the velocity potential and consequent section forces and moments 
for the aileron have been provided to the fifth power of a frequency 
parameter. In addition total force and moment coefficients have been 
obtained. 

The section coefficients were applied to a particular flutter study 
which attempted to compare: 

(a) The flutter characteristics for the aileron in various spanwise 

locations 

(b) The flutter characteristics of mass-balanced and mass -unbalanced 

ailerons 

(c) The flutter properties as determined from a three-degree-of- 

freedom problem with properties obtained from the two-degree- 

of-freedom subcases 

Because of the rather limi ted nature of this study (for example, no Mach 
number investigation was made, only two mass-ba l a n ce conditions were 
employed, and only one density ratio was employed) , the results of this 
investigation were not intended to be construed as conclusive but should 
aid in suggesting topics for future investigation. 



32 


NACA TN 3611-4 


With regard to the spanwise location of the aileron upon the flutter 
characteristics of the ving, no conclusion can readily be drawn other 
than that the inboard location of the aileron is best from a flutter 
standpoint. The possible detrimental effect of applying mass balance 
to the control surface in this Mach number region, shown in NACA TN 3160, 
was also indicated in this study, especially where three degrees of 
freedom were considered. Where finite- span considerations were applied 
to the two-degree -of- freedom subcases, mass balance appeared to be ben- 
eficial from a flutter standpoint for the wing-bending — aileron case 
but not for the wing-torsion — aileron case. The results of the two- 
degree-of- freedom subcases appeared to reflect the flutter properties 
obtained from the three -degree-of -freedom problem. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 9 , 19 ^ 6 . 
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APPENDIX A 

COEFFICIENTS ASSOCIATED WITH WING MOTION 


Most of the egressions in this appendix have been taken from ref- 
erence b, where they were presented as a power series containing terms 
up to the seventh power of the frequency parameter o>. Inasmuch as the 
present paper considers expansions to the fifth power of the frequency 
parameter, the expressions of reference ^ are accordingly abridged. These 
abridged expressions are repeated here for the sake of completeness. 


Section Coefficients Associated With Wing Motion 

For sections between the wing tip and the point where the leading- 
edge Mach line intersects the trailing edge, 0 < y < l/p, the section 
coefficients are: 




F i - F 2 


1 + - 6^9, + - x)f J 

r 3e L j 




l 2 „ A. feb + ^|(2P 2 - l)f 2 - *^ 3 ] + ^—[^2 - to* - ♦ 5)\ + 
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A + x[(6p^ + #2 - 1)9 2 - **to 2 + 5)F 5 ] + ^[(P a + - 

ac 2 1 p *L ^ 6p 8 L 
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In equations (A2) 


*n - Fn(l,y) 


For sections inboard of the point where the leading-edge Mach line 
intersects the trailing edge, l/f3 < y, the section coefficients are: 
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J pk 2 3p5 6op 9 


lx „ _A(s 2 - 1) + !^(P 2 + 5) - + >*9P 2 + 63 ) - 2xoL2 = 14’ - 2XOL2 1 

p3k 6p7 lSopJ-i 
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*iM 2 k 2 / 
I5p7 

= M-l' - 2x 0 L 1 ' 


% = i? - + 


m 2 =. A - 2£>e + 2*0 /,| C ? 

pk ap5 + 36p9 W + 7) - ^ 


= M 2 ‘ - 2XOL2 1 


“ 3 = A - i (p2 + 31 + + 3 =® 2 + 35) - a^Oij. + Lj + a^) 


pk 2 2p5 j6p9 

= M 3 * - axoCMi* + L 5 ' - 2x0^*) 


Mk = ^ (f2 ' X) + + 5) - ~rf^ k + *^ 2 + 63 ) - 


3p^k 15P 7 ' " 105P 11 ' 

2xo(M 2 + Llj. + 2 xoL 2 ) 

= Mj,.' - 2x 0 (M 2 I + L 4 ' - 2xoL 2 ') 


Total Coefficients Associated With Wing Motion 
^ >oTo!\r p “ ' = 2 *“ ° f ^3) not pro^ in 


> w 
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ia 


l 2 


3 - M 2 ^ 2 .( 4 p 2 + 5 ) 


p 3 6 p 7 


_l_ 

3A 


gl+2 _ tf^ p 4 + 8p 2 + 8) I 


r 5 P C 


+ M ^^ 2 + 7) 

P k p5 ipp9 


+ 12 p 2 + l£) 


_1_ 

2A 


_i_ - M!|(g 2 + 4 ) + 
p 2 * 3 ?® 


45 P 


10 


-i i-(p 2 + 3) + ^(4^ + 35P 2 + 35) 


[pk 2 3P 5 


6op 9 


_ 1 _ 

2A 


p 2 ^ 2 


-i-<P 2 + 4 ) + 3 ^^( 5 P 1< ' + 20 P 2 + 16 ) 
3P 6 45g 10 


- 


l 4 


-Up 2 . x) + J£W + 5 ) . (4p4 + 49P 2 + 63) 

7 iSop 11 


3 

pTt 


6 p 


3A 


p^ht 


ls^(p2 * 2) + + l6p 2 + 16) 


5 P l 


M 5 p- 


12 


- SXoLa 


(A5) 


Mn 


m 2 = 


it- -i<^4p 2 + 5) 

i 3 15P 7 


3P 


1 _ 

3A 


I 5 (p 2 + 2) . «^<p4 + 8p 2 + 8) 
[ 2 ^ 3P° 


- 2 XqLi 


i - + 2M(4p2 + 7) 

0 k 2 p 5 36 p 9 


_ 1 _ 

3A 


— - ^(p 2 + 4) + 
P^k 5P 6 


W^p4 + i2p 2 + 16) 
35P W 


- 2 XoL 2 


Mst 


U Up 2 + 3) + ?iM(4p4 + 35P 2 + 35) 

I pit 2 sp5 36p9 


5A 


2 

P 2 ^ 


JL_(3P 2 + 4) + i&S^(5P 4 + 20p 2 + l6 ) - a*b(% + c 3 + 2 'o&l> 
5P 6 35P W 


1 % 


5 P 

JL-(p 2 - i) + ^p 2 + 5) - 44p^ + 4 9 p 2 + 63) 

3p 3 k 15P 7 105P 11 — 

- 2 x 0 ( 55 2 + I4 + 2X0L2) 


+ i 

3A 


3 

p 4 k 


4M^( p2 + 2) + (3p4 + l6p 2 + l6) 

3P 8 15P 12 
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(a) Coefficients contributing to wing lift force and moment. 

* 

Figure 3 .- Spanwise distribution of components of lift-force and moment 
coefficients introduced by the presence of a control surface. M = 1.3 
k = l/9j A = 4j and y-j_/b = 2. 
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(a) Aileron at inboard location. 



(b) 


Aileron at mid span location. 



(c) Aileron at outboard location. 

Figure 4.- Flutter parameter b%/c plotted against frequency ratio 
0 ) 5 / 1 % for various aileron locations. Three degrees of freedom; 

M = 1.3j x p - 0.2} r^2 = 0.125. (Other parameters for wing and 
aileron are given in the text. ) 
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Calculated values 

Wing -bending— torsion values 

Estimated values 
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(a) Aileron at midspan location. 


Unstable 
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(To) Aileron at outboard location. 

Figure 5. - Flutter parameter bc%/c plotted against frequency ratio 
for various aileron locations. Three degrees of freedom; 

M = 1.3; xp = 0; r ^ = 0.100. (other parameters for wing and aileron 
are given in the text. ) 
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(a) Aileron at inboard location. 



(b) Aileron at midspan location. 
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(c) Aileron at outboard location. 

Figure 6.- Flutter parameter bo^/c plotted against frequency ratio 
cug/co^ for wing -bending -—aileron case employing finite-span theory. 

M = 1.3- 
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(a) Aileron at midspan location. 


Unstable 


(b) Aileron at outboard location. 

Figure 7«“ Flutter parameter botc/c plotted against frequency ratio 
ag/cjOh for wing- torsion— aileron rotation case employing finite-span 

theory. M = 1.3? — = 2.08 S2& = 0.480 
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